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It is now an accepted fact that the size at which dunes form from a fiat sand bed as well as 
their 'minimal size' scales on the flux saturation length. This length is by definition the relaxation 
length of the slowest mode toward equilibrium transport. The model presented by Parteli, Duran 
and Herrmann [Phys. Rev. E 75, 011301 (2007)] predicts that the saturation length decreases to 
zero as the inverse of the wind shear stress far from the threshold. We first show that their model 
is not self-consistent: even under large wind, the relaxation rate is limited by grain inertia and thus 
can not decrease to zero. A key argument presented by these authors comes from the discussion of 
the typical dune wavelength on Mars (650 m) on the basis of which they refute the scaling of the 
dune size with the drag length evidenced by Claudin and Andreotti [Earth Pla. Sci. Lett. 252, 30 
(2006)]. They instead propose that Martian dunes, composed of large grains (500 fim), were formed 
in the past under very strong winds. We show that this saltating grain size, estimated from thermal 
diffusion measurements, is not reliable. Moreover, the microscopic photographs taken by the rovers 
on Martian aeolian bedforms show a grain size of 87 ± 25 fim together with hematite spherules at 
millimetre scale. As those so-called "blueberries" can not be entrained by reasonable winds, we 
conclude that the saltating grains on Mars are the small ones, which gives a second strong argument 
against the model of Parteli et al. . 

PACS numbers: 45.70.Mg, 83.50.Ax 



In the present comment, we adopt the point of view 
of Parteh, Duran and Herrmann P] and use their model 
to point inconsistencies. We refer the interested reader 
to the series of papers published by the authors on the 

subject miaiiiisiEiiTiiH]. 



MODELLING THE SATURATION LENGTH 

The sand transport model used in [1 belongs to the 
series of models - the 'one species' models - in which one 
assumes that there is a single type of grain trajectories. 
The only self-consistent model of this type [4J is that 
derived by Ungar and HafF [9], from which [1] is directly 
inspired. One assumes that the evolution of the sand 
flux is governed by the ejection process. Introducing the 
grain hop length i and the number of ejected grains N 
per unit impacting grain, one gets: 



dx 



Nq. 



(1) 



The fluid in the saltation curtain is assumed to be at 
equilibrium between the driving shear stress PfU^, the 
air-borne basal shear stress PfU^^^ and the sand-borne 
shear stress. The sand-borne shear stress is proportional 
to the sand flux and to the difference between the velocity 
at which grains take off and collide back the sand bed 



At saturation, the wind is assumed to be just sufficient to 
maintain transport {N = 0), which leads to a basal shear 
velocity Wbas independent of and thus equal to the 
threshold shear velocity uth- At saturation U| and W| 
are evaluated in the saltation curtain, where the velocity 
proflle is almost independent on u*. The saturated flux 
can thus be put under the form: 



(3) 



where x depends only on the grain size, for a given at- 
mosphere, but not on the wind strength. 

Parteli et al. then derive the saturation length by 
a simple linearization of the saturation equation under 
three assumptions. First, the number of ejected grains 
N per unit impacting grain is assumed to be a function of 
the basal shear velocity Ubas only. Second, they assume 
that the grains ejected during collisions instantaneously 
reach the wind speed in the saltation curtain: the grains 
are assumed to have negligible inertia. Third, they as- 
sume that the wind instantaneously adjusts to changes 
of sand flux. One then gets: 
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This equation can be put under the form of a flrst order 
linear relaxation: 



' dx 



Qsat - q, 



(5) 



where the saturation length is equal to: 
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FIG. 1: Trajectory of a grain of 250 fim in the saltation cur- 
tain, on Earth. The transient length allows to define and 
measure the drag length fdrag- 



Parteli et al. have estimated the prefactor (-/ for a 
grain size of 250 /im, on Earth, to 0.85 m. As expected 
for any relaxation length, ^sat diverges at the threshold 
shear velocity. As £ does not depend on , tgat decreases 
as l/u1 for large u*. 

In reality, there is not a single mechanism limiting the 
time and length of saturation transient but several : (i) 
the ejection of grains, (ii) the grain inertia that controls 
the length needed for one grain to reach its asymptotic 
trajectory, (iii) the fluid inertia that controls the length 
needed for the wind to readapt to a change of q, (iv) the 
presence of grains above the saltation curtain with much 
longer trajectories It is worth emphasising that the 
ejection of grains (i) is the single source of lag considered 
in [T]. One should consider for 4at, the slowest process 
i.e. the largest relaxation length amongst the modes of 
relaxation. We first wish to show that the saturation 
length proposed by Parteli et al. is smaller than the re- 
laxation length imposed by the grain inertia. The equa- 
tion governing the grain motion may be written under 
the form: 




where the drag coefficient is approximated by 

2 

. Coo is the drag coefficient in the fully developed turbu- 
lent regime, i.e. at large particle Reynolds number. In 
this limit, the drag length ^drag, defined as the length 
needed for the grain to reach its asymptotic velocity, 
scales as ps/p/d. Consistently with Parteli et al., one 
can use for natural sand grains [5J [TU] Coo — 1 and 
s ~ 5. Reasonable collisions rules are those considered 
in [B in H], with a restitution coefficient e ~ 0.6 and a 
rebound angle around 45° . Figure [l] presents the tra- 
jectory of an ejected grain submitted to a wind at Uthi 
together with a fit of its envelop by an exponential relax- 
ation. The grain size is chosen to d = 250 pm, for the 



seek of comparison with [I]. We find a drag length, of the 
order of 570 mm. This means that the relaxation length 
associated to the ejection mechanism becomes smaller 
than that governed by grain inertia at moderately large 
velocity (around m, = 1.5Mth in fig- [2]). 

The model of Parteli et al. can be slightly modified to 
introduce the lag between the ejection of grains and the 
point at which they reach the saltation curtain velocity: 



dx^ and 4rag-^ 



(8) 



where Q is the flux of grains just ejected and already 
accelerated by the wind. The saturation length, defined 
for this second order system as the slowest relaxation 
rate, then becomes: 



dN 



A) 



(9) 



It is plotted in figure [2] together with the prediction by 
Parteli et al.. One can see that the divergence of the 
saturation length at the threshold is due to the ejection 
process, as stated by Parteli et al.. However, soon above 
the threshold (above — l.lSuth in fig- [2]), there is at 
least another mechanism leading to a larger saturation 
length: the grain inertia. 

We reach the first conclusion of this comment: as 
the saturation length cannot be smaller than the drag 
length, it cannot decrease with the wind strength far 
from the threshold. This is an evidence for the lack of 
self-consistence of the model proposed by Parteli et al.. 
The grain inertia could well be the limiting mechanism 
at large wind, as proposed in [H |3J [TT], but this does 
not preclude the existence of even slower relaxation pro- 
cesses mil]. 



THE SIZE AND DENSITY OF GRAINS ON 
MARS 

The main argument presented by Parteli et al. in 
favour of a saturation length decreasing as the inverse 
of the wind shear stress comes from the typical size of 
Martian dunes. They use an estimated grain size d ~ 
500±100 pm derived by Edgett and Christensen [T2 from 
the Viking Orbiter infrared thermal mapper (IRTM) 
data. Using the simple scaling law of the dune wave- 
length based on .^drag [H]) one would then expect with 
such a grain size a spacing of 4 km between dunes on 
Mars. The real wavelength of Martian dunes is much 
smaller, between 500 m and 700 m |8]. Parteli et al. thus 
conclude that, in disagreement with our scaling relation- 
ship, very large winds are needed to explain the observed 
sizes i.e. to make the saturation length of large grains 
very small. We hereafter summarise the evidences given 
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FIG. 2: Saturation length as a function of the rescaled wind 
shear velocity for sand grains of 250 fim, on Earth. The dotted 
line corresponds to the model of Parteli et al. (the symbols 
are deduced from fig. 5 of [l]), that only takes into account 
the lag due to the ejection of grains. The solid line is the 
relaxation length obtained by modifying the model to take 
into account the grain inertia (Eq. [8| . The sharp transition is 
to be related to the second order dynamics. The dashed line 
shows the value of the drag length. 



in [8] that the grains in saltation on Mars are in fact 
much smaller than 500 fim. 

As shown by Fenton et al. [ISlIT^, the determination 
of saltating grain size from thermal diffusion estimates 
is far to be obvious. Indeed, the measurement is very 
indirect. Edgett and Christensen \IT use the thermal 
model of Kieffer et al. [15j to calculate thermal inertia 
with Viking IRTM data. Using the updated relation by 
Presley and Christensen [TB^ , the same data for the HcUe- 
spontus dunes give 1200±200 fim instead of 500±100 /im. 
If we push further the argument of Parteli et al. with this 
new grain size, one would need a typical shear velocity 
of 10 m/s, which roughly corresponds to 500 km/hour 
at 10 m above the soil. This is more than one order of 
magnitude larger than present winds observed on Mars. 
In a dust devil rotating at such a speed, the depression in 
the core of the vortex would be larger than the average 
pressure of Martian atmosphere, which is not physically 
possible. 

Much more reliable are the direct observations by the 
rovers Opportunity and Spirit. The photographs taken 
by the rovers (Figure |3^-b) mostly show two well sep- 
arated grain sizes: large spherules of millimetric scale, 
composed of hematite {ps — 5270 kg/m'^) and small 
basalt grains with iron coating {ps = 3010 kg/m'^) be- 
tween 60 iim and 110 /xm [S1[T7]. Note that the thermal 
inertia measurements would point to large grains even in 
the zones where the rovers have found small ones [18]. 
How can one then discriminate between grains that can 
be transported in saltation and grains that cannot? The 
first argument is theoretical [8 . With such large size and 



density, the threshold velocity for the entrainment of the 
blueberries into saltation is very large (fig.|4]). The small 
grains, on the other hand, can be transported even with 
contemporary winds. It may be argued that the small 
grains would rather be transported into suspension but 
there is no clear threshold between saltation and suspen- 
sion: as the wind speed increases, wind fluctuations be- 
come gradually more important with respect to gravity. 
The second argument comes from the field observations of 
aeolian structures on Earth. In particular, the formation 
of aeolian ripples and of shadow dunes behind obstacles 
(nebkhas) constitute a clear signature of transport into 
saltation. Figure [3j:-e shows that these structures have 
a much higher concentration of small particles than the 
surrounding soil. 

The emergent picture is thus very coherent: the grains 
transported in saltation on Mars are smaller than 100 /im 
and are certainly not the millimetre scale hematite 
spherules blueberries; they can be transported by the 
present winds (ripples have formed in very recent impact 
craters; they form very recognisable aeolian bedforms like 
ripples and nebkhas, and probably dunes. The conclu- 
sions reached by Parteli et al. are thus probably wrong, 
for the particular problem of dune formation on Mars 
and for the modelling of sand flux saturation transients 
in general. 



RELATION BETWEEN THE WAVELENGTH AT 
WHICH DUNES FORM AND THE SATURATION 
LENGTH 

The instability of a flat sand bed results from the inter- 
action between the sand bed profile, which modifies the 
fluid velocity field, and the fiow that modifies in turn the 
sand bed as it transport grains. The fiuid is accelerated 
on the upwind (stoss) side of proto-dunes and deceler- 
ated on the downwind side. This results into an increase 
of the shear velocity applied by the fiow on the stoss 
side of the bump. Conversely, decreases on the lee 
side. As the saturated sand fiux is an increasing function 
of u*, erosion takes place on the stoss slope as the fiux 
increases, and sand is deposited on the lee of the bump. 
If the velocity field was symmetric around the bump, 
the transition between erosion and deposition would be 
exactly at the crest, and this would lead to a pure prop- 
agation of the bump, without any change in amplitude 
(the 'A' effect). In fact, due to the simultaneous effects of 
inertia and dissipation, the velocity field is asymmetric 
(even on a symmetrical bump) and the position of the 
maximum shear stress is shifted upwind the crest of the 
bump (the 'B' effect). In addition, the sand transport 
reaches its saturated value with a spatial lag ^sat- The 
maximum of the sand fiux q is thus shifted downwind 
the point at which is maximum by a typical distance 
of the order of ^sat- The criterion of instability is then 
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FIG. 3: a Microscope photograph of the sand composing a Martian ripple, b Microscope photograph showing the mixing of 
small grains and hematite spherules ("blueberries"), characteristic from the soil seen by the two rovers, c Aeolian ripple on 
Mars, characteristic of transport in saltation. A strong difference of composition between the soil, covered by blueberries and 
the ripple can be observed, d Aeolian nebkhas on Mars, characteristic of transport in saltation. These shadow dunes behind 
stones are clearly evidencing that small grains are transported in saltation, but not the hematite blueberries, e Extended 
zone of aeolian ripples in a small scale impact crater. Blueberries may be seen at the bottom left of the picture, showing that 
the ripples are composed of small grains. These pictures have been taken by the rover Opportunity on Sols 58, 59, 60 and 
85. They can be found on the NASA web site http://marsrovers.jpl.nasa.gov/gallery/all/opportunity.html. Courtesy 
NASA/JPL-Caltech. 



geometrically related to the position at which the flux is 
maximum with respect to the top of the bump: an up- 
shifted position leads to a deposition of grains before the 
crest, so that the bump grows. 

These arguments can be formalized by performing the 



linear stability analysis of a flat sand bed [5] . For a 
small deformation of the bed profile h{t,x), the excess 
of stress induced by a non-flat profile can be written in 
Fourier space as pful{A + iB)kh. A and B may be in 
principle deduced from a turbulent closure. Jackson and 
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FIG. 4: Diagram showing the mode of transport on Earth (a) 
and on Mars (b) , as a function of the grain diameter d and of 
the turbulent shear velocity Uth- Below the dynamical thresh- 
old (dashed line), no grain motion is observed (dark gray). A 
grain at rest on the surface of the bed starts moving, dragged 
by the wind, when the velocity is above the static threshold 
(solid line). Between the dynamical and static thresholds, 
there is a zone of hysteresis where transport can sustain due 
to collision induced ejections. Above static threshold, the 
background color codes for the ratio u,/ufaii: white corre- 
sponds to negligible fluctuations and gray to suspension. The 
experimental points are taken from (O) Chepil [T5] and (□) 
Rasmussen [H |20] in the aeolian case. The insets shows the 
location of the observed grains in the diagrams. 



Hunt [2TJ 122] have derived asymptotic expressions for A 
and B as functions of In(fczo), where zg is the aerody- 
namic roughness. 

At this stage, for who wishes to catch subdominant 
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FIG. 5: a Relation between the wavelength Am at which dunes 
appear and the saturation wavelength, as a function of the 
wind strength. The limit of stability Ac (wavelength for which 
CT = 0) is also shown. The plot has been produced for A = 5 
and B = 1.5, which are typical values predicted by Jackson 
and Hunt |21II22| . The slope effect, due to gravity, has an in- 
creasing importance close to the threshold, b Destabilization 
of the flanks of a barchan dune during a violent dust storm 
coming from Sahara towards Canarias, in April 2003. The 
wavelength of destabilization is reduced by a rough factor of 
2 with respect to that observed during regular trade winds 
(~ 20 m). 



dependencies on the wind shear velocity, there is again 
a very important mechanism forgotten in Parteh et 
ai: the influence of the slope on the threshold shear 
stress. As shown by Rasmussen et al. [5U], at the lin- 
ear order, the threshold shear stress may be written as: 
PfU^Yi{\ + dxh/ tan 6'a), where Oa — 32° is the avalanche 
repose angle. This gravity effect originates from the trap- 
ping of grains at the surface of the sand bed, influenced 
by the slope. In the Fourier space, the saturated flux 
modulation can be written as: 



9sat = [{A + iB)ul 



lIJ tan( 



Xkh. (10) 



Using a first order linear saturation equation and the 
conservation of matter, we end with a growth rate a of 
the form: 
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Figure [s] shows the relation between the wavelength Am 
at maximum growth rate and the marginally stable wave- 
length Ac, as functions of the rescaled wind shear velocity 
w*/uth. It can be observed that these wavelengths de- 
crease with wind strength, due to the decreasing relative 
importance of gravity effects with respect to wind effects. 



This slope effect could be in fact the dominant explana- 
tion for the observed variations of minimal size with the 
wind strength (typically a factor of 2 in Morocco, see 
figure [5]). 

In conclusion, if the role of the particle diameter and 
of the fluid to grain density ratio on the time and length 
scales of dunes is now pretty clear, that of wind speed re- 
mains controversial. Further work is needed to shed light 
on the influence of the numerous dynamical mechanisms 
involved. 
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